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We present precision calculations of the processes e + e~ — > 4-fermions 
in which the double resonant W + W~ and ZZ intermediate states occur. 
Referring to these latter intermediate states as the 'signal processes', we 
show that, by using the YFS Monte Carlo event generators YFSWW3-1.14 
and KoralW-1.42 in an appropriate combination, we achieve a physical 
precision on the WW signal process, as isolated with LEP2 MC Workshop 
cuts, below 0.5%. We stress the full gauge invariance of our calculations 
and we compare our results with those of other authors where appropri- 
ate. In particular, sample Monte Carlo data are explicitly illustrated and 
compared with the results of the program RacoonWW of Denner et al. . In 
this way, we cross check that the total (physical©technical) precision tag 
for the WW signal process cross section is 0.4% for 200 GeV, for example. 
Results are also given for 500 GeV with an eye toward the LC. For the 
analogous ZZ case, we cross check that our YFSZZ calculation yields a 
total precision tag of 2%, when it is compared to the results of ZZTO and 
GENTLE of Passarino and Bardin et al, respectively. 
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1 Introduction 



The theoretical paradigm affirmed by the award of the 1999 Nobel Prize to G. 't 
Hooft and M. Veltman for the success of the predictions of their formulation [Q] of the 
renormalised non-Abelian quantum loop corrections for the Standard Model [[| of the 
electroweak interaction focuses our efforts on the need to continue to test this theory 
at the quantum loop level in the gauge boson sector itself. This then emphasises 
the importance of the on-going (the data are under analysis and will be for some 
time even though the LEP2 accelerator was recently shutdown) precision studies of 
the processes e + e~ — > W + W~(ZZ) + 72(7) — ► 4/ + 71(7) at LEP2 energies 0,g,§, as 
well as the importance of the planned future higher energy studies of such processes 
in LC physics programs |6|,|7|,§,[| . We need to stress also that hadron colliders also 
have considerable reach into this physics and we hope to come back to their roles 
elsewhere pD|. 



In what follows, we present precision predictions for the event selections (ES) of 
the LEP2 MC Workshop [[TIJ for the processes e + e~ -> W + W~ + 72(7) -> 4/ + 71(7) 
based on our new exact 0(a) pro d YFS exponentiated LL 0(a 2 ) FSR leading pole 
approximation (LPA) formulation as it is realized in the MC program YFSWW3- 
1.14 |jl2|,|T3l in combination with the all four-fermion processes MC event generator 



KoralW-1.42 | | so that the respective four-fermion background processes are taken 
into account in a gauge invariant way. In addition, we also present the current status 
of the predictions of our YFS MC approach to the processes e + e~ — > ZZ + 72(7) — > 
4/ + 77(7) as it was also illustrated in the 2000 LEP2 MC Workshop [ ] using the 
MC event generator YFSZZ @, which realizes YFS exponentiated LL 0{a 2 ) ISR 
in the LPA in a gauge invariant way. Indeed, gauge invariance is a crucial aspect 
of our work and we stress that we maintain it through-out our calculations. Here, 
ISR denotes initial state radiation, FSR denotes final-state radiation and LL denotes 
leading-log as usual. 

This realization which we present of the YFS MC approach is the exclusive expo- 
nentiation (EEX) |L6| and it is already well established in its applications to the MC 
event generators for LEP1 physics calculations in the MC's KORALZ/YFS3 |17,[18], 
BHLUMI [|Tggg] and KoralW 0. In our applications in YFS WW3- 1.14 and in 
KoralW-1.42, the FSR is implemented using the program PHOTOS so that not 
only is the FSR calculated to the LL 0(a 2 ) but the FSR photons have the correct 
finite pt in the soft limit to 0{a). We always use the ratio of branching ratios (BR's) 
to correct the respective decay rates through O(a) accordingly. Recently, we have 
introduced the coherent exclusive exponentiation (CEEX) approach to the YFS 
MC event generator calculation of radiative corrections and we will present the ap- 
plication of this new approach to the 4/ production processes elsewhere |T0]]. For a 
description of its application to the 2f production processes see Ref. [23 . 

Recently, the authors in Refs. [p§]25[ have also presented MC program results for 
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the processes e + e~ — *> W + W~ + 71(7) — > 4/ + 71(7), n = 0, 1 in combination with 
the complete background processes which feature the exact LPA 0{a) correction, the 
complete 0(a) result for e + e~ — ■> 4/ + 7, and soft photon KF | § exponentiation for 
the LL 0(a 3 ) ISR via structure functions. Thus, we will compare our results where 



possible with those in Refs. |24| in an effort to check the over-all precision of our 
work. As we argue below, the two sets of results should agree at a level below 0.5% 
on observables such as the total cross section. The authors in Refs. have used semi- 
analytical methods to compute the exact LPA 0(a) correction e + e~ — > W + W~ — > 4/ 
with no higher order resummation. Thus, while we have compared our results with 



theirs in Ref. [13] for example, here we do not present such comparisons because the 



expected precision tag of their results is larger than the desired 0.5% needed by the 
LEP2 experiments [[Tl"|j . 

For the processes e + e~ — > ZZ + (wy) — > 4/ + (n'j), the authors in Refs. p7| . [28f 
have presented calculations in the LEP2 MC Workshop |lTJ at the NC02 and all-4/ 
level I ] as well. The calculations in Ref. | |] are done with the program ZZTO 
and feature universal ISR corrections, 0(a) FSRqed corrections, 0(a s ) FSRqcd 
corrections, and running masses in the fermion loop scheme of Ref. . The results 
in Ref. | | feature the structure function approach to the ISR QED corrections and 
the 0(a) FSRqed corrections. We will compare our YFSZZ results with these two 
sets of results as well, as the three approaches should agree at the level of the 2% 
precision needed by the LEP2 experiments [[II] on observables such as the total cross 
section. 

Our presentation is organised as follows. In the next Section, we discuss the cur- 
rent status YFSWW3-1.14. In Section 3, we present the current status of KoralW-1.42 
from the standpoint of its use to calculate the 4/ background processes in combina- 
tion with YFSWW3-1.14. In Section 4, we present the current status of YFSZZ. In 
Sections 5 , 6 and 7, we illustrate the results we have obtained with our calcula- 
tions for YFSWW3, KoralW-1.42 and YFSZZ, respectively, for the ES of the LEP2 
MC Workshop [[□]], wherein we include comparisons with the respective results in 
Refs. p4] , |27| , p8fl . Section 8 contains our summary remarks. 
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2 YFSWW3-1.14 



In this section we present the current status of YFSWW3-1.14. We start with the 
process of interest and its cross section, 

e-fa) + e + (p 2 ) - A(n) +7 2 (r 2 ) + f[(r[) + f 2 (r' 2 ) + 7 (fc : ), ..., 7 (M, 

o"n = 77— / ^T n+4 (pi + P2; ri, r 2 , r' l5 7-3, fci, fc„) /, v 

flux J {1) 

Y Y l-^4/(Pi»P2> r i,^2,^,r2,/ci,...,A; n )| 2 , 

ferm. spin phot, spin 

and the corresponding expressions for the W + W~ production and decay in 
the leading pole approximation (LPA), 

e"(pi) + e + (p 2 ) -> H/-( gi ) + W + (g 2 ), 

- A(rO +7a(r 2 ), ^ + fe) - AW) + 7' 2 (r 2 ), 

dr n+4 {pi +p 2 ;r 1 ,r 2 ,r' 1 ,r' 2 , h,...,k n ) (2) 



<7, 



flux 



r li r 2i ^1) ■•■) ^re) I 

ferm. spin phot, spin 



Here, we realize the LPA afc as follows: 



M < f f \pi,p 2 ,r 1 ,r 2 ,r' 1 ,r' 2 ,ki,...,k n ) L ^ M i [ l jl A (pi,p 2 ,r 1 ,r 2 ,r' 1 ,r 2 , k 1 , k n 

( (n),AiA 2 , 



Y M Prod 2 (pi,P2,Ql,q2,k 1 ,...,k a ) 

Phot. Partitions 

Djm) M Dlc 1 ,xM> ri > r2 > ka + 1 >--' kb ) (3) 



■ ,(n) 



D( qi ) = q 2 - M 2 , M 2 = {M w 2 - iT w M w ){l - T w 2 /M w 2 + 0(a 3 )), 
q l= n+r 2 + k a+ i + ... + k b ; q 2 = r' l +r' 2 + k b+1 + ... + k n , 

where the two formulations of the LPA, LPA a are based on the results in 
Eden Refs. ||30| , |31|| as one can see from the representation of our amplitudes 

M as 

M = Y J ^ 3 ({q k q l }). (4) 

3 

Here, the {£j} are a complete set of spinor covariants and the {A,} are 
the respective scalar functions. For LPA( a ) b , we do (not) evaluate the 
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spinor covariants on-pole in realizing the respective Ai^ PA . We do both in 
YFSWW3-1.14. 

We use standard YFS methods (EEX-Type) to write 



da 



maB'+2aB 



1 



{2-kY 

d 3 r± d?r 2 d 3 r[ d 3 r' 2 
X E 1 E 2 E[ % 



n=l i 



j Pn{kl, k n )] 



(5) 



where 



D 
2aB 



d 3 k^ 



fk 

kn 



S 



-iyk 



e(K n 



\k\) 



(6) 



9{K max - \k\)S{k). 



Here, K max is a dummy parameter of which eq.@ is independent. In 
realizing eq.(|5|) in YFSWW3, we employ the following schemes, which are 
related by the renormalization group: 



• Version 1.13: G^-Scheme of Fleischer et al. 

• Version 1.14: Scheme A — only the hard EW correction has «g m ; 
Scheme B — the entire O(a) correction has a(0) 

As it was shown in the LEP2 MC Workshop [|11[], there is a -0.3 -=--0.4% 
shift of the normalisation of version 1.14 relative to that of version 1.13. 
This can be seen as follows. The universal LL ISR O(a) soft plus virtual 
correction is 



a 



7T 



7T 
~3 



(7) 



with P=^(L — 1) and with kn equal to the usual soft cut-off and L = \\is/m 2 e . From 
eq.(^), we get the estimate of the shift in normalisation between version 1.13 and 
version 1.14 at 200 GeV as 



(a(0)-oo M )(-L-2)~-0. 



1 . 



This is consistent with what is observed as reported in Ref. |Ll]. See Dittmaier's 
talk p5i for more details and references. 
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3 KoralW-1.42 



For the process of interest, e~(p 1 ) + e + (p 2 ) — > + f 2 {r2) + f[( r i) + f 2 { r ' 2 ) + 

7(^1), j(k n ), we use KoralW-1.42 which realizes the 0(a 3 ) LL YFS expo- 
nentiated ISR. The respective input Born matrix elements are the GRACE 



v. 2 [|33|| all 4f library of Born matrix elements and our independent CC03 
Born matrix elements. This allows us to combine YFSWW3-1.14 and 
KoralW-1.42 to correct for background diagram effects: using LPA„ in 
YFSWW3-1.14, whose cross section we denote by a(Y a ), we get 

a Y/K = a(Y a ) + Aa(K), (9) 

where Aa(K) is defined by 

Aa(K) = a(K 1 )-a(K 3 ). (10) 

Here, a{K x ) is the 4-f KoralW-1.42 result and a(K 3 ) is the CC03 KoralW- 
1.42 result. This means that ay/K 1S accurate to 0{^j^-). 

Alternatively, using LPA i5 i — a,b in YFSWW3-1.14, whose cross sec- 
tion we denote by o~(Yj), we get 

a K/Y = a(K l ) + Aa(Y) (11) 

where 

Aa(Y) = a(Y i )-a(Y 4 ), (12) 

and 0(^4) is the respective YFSWW3-1.14 result with NL O(a) corrections 
to P n , n — 0, 1, switched off. This means that <Jkiy is also accurate to 

Above WW threshold, okiy and ay/K 
agree to the 0.1% level. We advocate the latter as our best result in the 
following. 

Note that we sometimes identify cr(Yi) = cr(Y a ), a(Y 2 ) = 0"(F&), cr(Y 3 ) = 
cr(K 3 ) with a(A' 2 ) equal to the cross section from KoralW-1.42 with the 
on-pole CC03 Born level matrix element with YFS exponentiated 0(a 3 ) 
LL ISR — this a(K 2 ) should be available soon. It is useful for further cross 
checks on our work. 



4 YFSZZ 

In our calculation in YFSZZ-1.02 |15| the process of interest is e~(pi) + e + (p 2 ) 
-> Z( qi )Z(q 2 ) + (7^),. -,7(U) -/i(ri) +/i(r 2 ) + ft(r[) +j' 1 (r' 2 )+ 1 (k 1 ),..., 1 (k n ). 
We proceed as follows in realizing the MC YFSZZ-1.02: 
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• We use LPA a as described above for the NC02 process to calculate 
0(a 2 ) LL YFS exponentiated ISR for the input NC02 Born matrix 
elements of Ref. [|34fl . 



• Anomalous couplings are supported following the conventions of Ref. | 34 | 
- this is also true for YFS WW3/KORALW. 

We stress that YFSZZ is in wide use at LEP and that it was tested in 
the LEP2 MC Workshop, just as YFSWW3-1.14 was tested. We now turn 
to such results. 



5 Results- YFS WW3- 1.14 



In this section we illustrate the effects of the NL O(a) correction as it is calculated 
in YFSWW3-1.14. We do this with the hardest photon angular distribution. Similar 
calculations of other observables can be found in Ref. 



ran 



Specifically, in Fig. p], we show the distribution of the cosine of the production 
angle of the hardest photon in the cms system with respect to the e + beam. We see 
that away from the beams the NL 0(a) correction is important for precision studies 
of this photonic observable. Similar conclusions follow from the more complete set of 
observables studied in Refs. [p^JTT| . 

Indeed, in the LEP2 MC Workshop, we compared our results with those of 
RacoonWW by the authors in Ref. p4| . For a complete description of these compar- 
isons we refer the reader to Ref. [11|. Here, we show in Table |l| the comparison of the 
total cross sections at 200 GeV with no cuts as defined Ref. iTTfl. 
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Hardest Photon Angular Distribution 
E CM = 200 GeV E CM = 500 GeV 



Aim 



• ISR: BARE 

* Best: BARE 
= ISR: CALO 
o Best: CALO 



-0.5 0.0 0.5 1 
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'i : : r.'. I '.^i . cos 

-0.5 0.0 0.5 1.0 



Figure 1: cos9-y w.r.t. the e + beam in the cms system for e + e — ► W + W — 
We see that NL corrections are important away from the beams, for example. 



ud(j, v 
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comparison with RacoonWW 



no cuts 


ctot [fb] 


final state 


program 


Born 


best 




YFSWW3 


219.770(23) 


199.995(62) 




RacoonWW 


219.836(40) 


199.551(46) 




(Y-R)/Y 


-0.03(2)% 


0.22(4)% 




YFSWW3 


659.64(07) 


622.71(19) 




RacoonWW 


659.51(12) 


621.06(14) 




(Y-R)/Y 


0.02(2)% 


0.27(4)% 




YFSWW3 


1978.18(21) 


1937.40(61) 


udsc 


RacoonWW 


1978.53(36) 


1932.20(44) 




(Y-R)/Y 


-0.02(2)% 


0.27(4)% 



Table 1: Total cross sections, CC03 from RacoonWW, YFSWW3, y/s = 200 GeV without cuts. 
Statistical errors correspond to the last digits in ( ). 



From the results in Table [I] and the related results given in Ref. [|IT|] we conclude 
that the TU of the calculations is 0.4% at 200 GeV for the total signal cross section. 
This is a considerable improvement over previously quoted precision of 2% in Ref. |jq] . 



6 Results- YFSWW3 /KoralW 

One of the important aspects of the isolation and study of the WW signal processes 
is the control of the corresponding background 4/ processes. This we do with our 
all-4/ MC KoralW-1.42 as we described above. Here, we illustrate the size of the 
corresponding 4/ background corrections to the YFSWW3-1.14 cross sections. 

Specifically, in Tabs. |2] and |3|, we show the size of this 4/ background correction 
in comparison to the NL correction of YFSWW3-1.14 for the total cross section, for 
example, both for the case of no cuts and the case of cuts, respectively, as defined in 



Ref. [ 1]. These results show that the 4/ background correction at 200 GeV to the 



total YFSWW3-1.14 cross section is below 0.1 ( 



7 Results- YFSZZ 

In this section we show the results of our comparison with ZZTO for the total 



ZZ pair signal processes as carried out in Ref. |y] . In that same set of comparisons 



ZZTO was also compared with the results of GENTLE by the authors in Refs. PS ]. 
In this way, a cross check was made on all three calculations. 
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WW/4f Cross Section 



NO CUTS 


&ww [fb] 


5 Af 


[%] 


°WW l /o l 


Final state 


Program 


Born 


ISR 


Born 






YFSWW3 


219.793(16) 


204.198 (09) 


— 


— 


-1.92 (4) 




KoralW 


219.766 (26) 


204.178 (21) 


0.041 


0.044 


— 




(Y-K)/Y 


0.01 (1)% 


0.01(1)% 










YFSWW3 


659.69 (5) 


635.81 (3) 






-1.99 (4) 




KoralW 


659.59 (8) 


635.69 (7) 


0.073 


0.073 






(Y-K)/Y 


0.02(1)% 


0.02(1)% 










YFSWW3 


1978.37(14) 


1978.00 (09) 






-2.06 (4) 


udsc 


KoralW 


1977.89 (25) 


1977.64 (21) 


0.060 


0.061 






(Y-K)/Y 


0.02(1)% 


0.02(1)% 









Table 2: Total WW YFSWW3 and KoralW cross sections: Born and ISR level, KoralW 
4/ correction, YFSWW3 0{a) NL correction, at 200 GeV , no cuts. The last digits in (• • • ) 
correspond to the statistical errors. 



Specifically, we show in Table 
predicted by YFSZZ and ZZTO for the case of no cuts as defined in Ref 



the ZZ signal cross section at 188.6 GeV as 

rj. For 



ZZTO, results are shown for two schemes, the and a schemes [^7|. The agreement 
between 

the programs in this comparison and between the programs in the other related 
comparisons carried out in Ref. [11] show that the TU for the respective NC02 signal 
process is 2% at the respective LEP2 energies. 



8 Conclusions 

We are currently at an exciting point in the tests of the EW Theory in gauge 
boson physics. The WW pair production is an important aspect of these tests. The 
radiative corrections which we realize in YFSWW3-1.14 play a significant role in these 
tests as follows: 

• Mass distributions: these are affected by FSR, yielding peak position and height 
shifts 

• W Angular distributions: these are affected by LL and NL corrections 

• £ Angular distributions: these are affected by LL and NL corrections 

• Photon Angular distributions: these are affected by LL and NL corrections 
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WW/4f Cross Section 



WITH CUTS 


&ww [fb] 


5 Af 


[%] 


°WW l /o l 


Final state 


Program 


Born 


ISR 


Born 






YFSWW3 


210.938(16) 


196.205 (09) 


— 


— 


-1.93 (4) 




KoralW 


210.911 (26) 


196.174(21) 


0.041 


0.044 


— 




(Y-K)/Y 


0.01 (1)% 


0.02(1)% 










YFSWW3 


627.22 (5) 


605.18(3) 






-2.00 (4) 




KoralW 


627.13(8) 


605.03 (7) 


0.074 


0.074 






(Y-K)/Y 


0.01(1)% 


0.02(1)% 










YFSWW3 


1863.60(15) 


1865.00 (09) 






-2.06 (4) 


udsc 


KoralW 


1863.07(25) 


1864.62 (21) 


0.065 


0.064 






(Y-K)/Y 


0.03 (2)% 


0.02(1)% 









Table 3: Total WW YFSWW3 and KoralW cross sections: Born and ISR level, KoralW 4/ 
correction, YFSWW3 0{a) NL correction, at 200 GeV, with cuts. The last digits in (• • • ) 
correspond to the statistical errors. 



Comparison with ZZTO 



channel 


YFSZZ 


ZZTO Gp -scheme 


ZZTO a-scheme 


qqqq 


294.6794(490) 


298.4411(60) 


294.5715(59) 


qqw 


175.4404(302) 


175.5622(35) 


174.9855(35) 


qq\\ 


88.1805(134) 


88.7146(18) 


87.9881(18) 


\\w 


26.2530(463) 


26.0940(5) 


26.1342(5) 


1111 


6.5983(15) 


6.5929(1) 


6.5706(1) 


WW 


26.1080(71) 


25.8192(5) 


25.9868(5) 


total 


617.2596(755) 


621.2241(124) 


616.2366(123) 



Table 4: NC02 cross sections, YFSZZ vs ZZTO, 188.6 GeV, in fb. The statistical errors 
correspond to the last digits in ( ). 

• Photon Energy distributions: these are affected by LL corrections 

• Normalisation: this is affected by LL AND NL corrections; 

the current 200 GeV TU is 0.4% from the {YFSWW3/RacoonWW} results. 

Concerning our results on calculating the 4/ background to YFSWW3-1.14 using 
KoralW-1.42, we have shown the following: 

• Two different combinations of YFSWW3 and KoralW-1.42 cross sec- 
tions reach the total precision C(^-^ SL ). 
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• The size of the 4/ correction to YFSWW3-1.14 is < 0.1%, as expected. 

• The future extension to a single platform is possible. 

It follows that YFSWW3/KoralW is a complete MC event generator 
solution for precision WW/4f production at LEP2 ( and LC's). 

From our studies of the NC02 signal process we conclude that YFSZZ, a multiple 
photon MC event generator for NC02 with j3 level LPA YFS exponentiation (EEX), 
is tested in the LEP2 MC Workshop vs ZZTO and GENTLE to 2% TU. An upgrade 
to higher precision is possible but is not needed, apparently? 
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